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Abstract 

The time evolution of excited nuclei, which are in equilibrium with the surrounding 
vapour, is investigated. It is shown that the finite nuclear systems undergo a first 
oder phase transition. The caloric curve is presented for excited 16 0, 24 Mg, 27 A1 



' and 40 Ca and the critical temperature is estimated for 16 0. 
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' 1 Introduction 



Mean-field models predict a first order phase transition for nuclear matter 



with a critical temperature which depends on the proton-neutron asymmetry 
[1—5]. A recent experimental attempt by the ALADIN group [6] to deduce an 
equation of state, which relates the excitation energy of a hot nucleus to its 
temperature, has stimulated both, theoretical and experimental efforts in this 
field [7-13]. In the experiment excited projectile spectators were investigated in 
Au+Au collisions at a beam energy of E/A = 600 MeV. While the equation 
of state refers to a stationary system where liquid and vapour (evaporated 
nucleons) are in equilibrium, the experiment deals with an expanding source. 
This causes some uncertainties for the temperature, which is deduced from 
isotope ratios, since the system cools while it is expanding. 

In molecular-dynamics calculations the finite system may be excited without 
flow, but similar problems arise when the phase transition sets in. Particles 
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which escape from the nucleus cool down the residue and thermal equilib- 
rium cannot be maintained. In order to avoid these difficulties, in the present 
simulations the excited nuclear system is confined by a wide container po- 
tential which is chosen to be a harmonic oscillator potential. Its frequency uj 
serves as a thermodynamic variable like the volume in the ideal gas case. Due 
to the containment evaporated nucleons cannot escape, but form a cloud of 
equilibrated vapour around the excited nucleus. 

Thermodynamic relations are obtained by coupling the nuclear system to a ref- 
erence system which serves as a thermometer. Both, the time-evolution of the 
nuclear system and of the thermometer are described by the Fermionic Molec- 
ular Dynamics (FMD) model. Assuming thermal equilibrium in the sense of 
ergodicity the temperature of the nuclear system is derived from the time- 
averaged energy of the thermometer and related to the excitation energy of 
the nucleus. 



2 Model and setup 

2. 1 The Fermionic Molecular Dynamics model 

The time evolution of the nuclear system is described within the framework 
of Fermionic Molecular Dynamics (FMD) published in detail in Ref. [14]. The 
model describes the many-body system with a parameterized antisymmetric 
many-body state \Q(t)) composed of single-particle Gaussian wave packets 



which are localized in phase space at f and p with a complex width a. Spin 
and isospin are chosen to be time-independent in theses calculations; they 
are represented by their z-components m s and m t , respectively. Given the 
Hamilton operator H the equations of motion for all parameters are derived 
from the time-dependent variational principle (operators are underlined with 
a tilde) 



( x | q(t) ) = exp 




<g> \m s )<S> | m t ) , 



(1) 
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In the present investigation the effective two-body nucleon-nucleon interac- 
tion V in the Hamilton operator consists of a short-range repulsive and long- 
range attractive central potential with spin and isospin admixtures and in- 
cludes the Coulomb potential [15]. The parameters of the interaction have 
been adjusted to minimize deviations between calculated and measured bind- 
ing energies and charge radii for nuclei with mass numbers 4 < A < 40. 

2.2 The container 

Besides the kinetic energy T and the nucleon-nucleon interaction V the Hamil- 
ton operator H includes an external field 



which serves as a container. 

The container is an important part of the setup because it keeps the evapo- 
rated nucleons (vapour) in the vicinity of the remaining liquid drop so that it 
equilibrates with the surrounding vapour. The vapour pressure is controlled 
by the external parameter u>, which appoints the accessible volume. 

2.3 The thermometer 

The concept of determining the temperature is to bring a reference system, for 
which thermodynamic relations between temperature and measurable quan- 
tities are known, into thermal equilibrium with the investigated system. The 
weakly interacting ideal gas, where the temperature is given by the mean ki- 
netic energy of the particles, may serve as an example. The reference system 
is called a heat bath if its heat capacity is much larger than that of the system 
and it is called a thermometer if its heat capacity is much less. 

As the nuclear system is quantal and strongly interacting its temperature can- 
not be deduced from the momentum distribution or the mean kinetic energy 
of the nucleons. Therefore, the concept of an external thermometer which is 
coupled to the nuclear system is used in the present investigation. The ther- 
mometer consists of a quantum system of distinguishable particles moving in 
a common harmonic oscillator potential different from the container potential. 

The time evolution of the whole system is described by the FMD equations of 
motion. For this purpose the many-body trial state is extended and contains 
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now both, the nucleonic degrees of freedom and the thermometer degrees of 
freedom 



\Q) = \Qn)® \Qth) ■ (4) 
The total Hamilton operator including the thermometer is given by 



H = H +H +H , H =T + V + U(uo) (5) 

with the nuclear Hamilton operator H and the thermometer Hamilton oper- 
ator 



/ ^ 2 
N th I 



H = E 

i=i 



2 m t/l («) 2 



fc (2) 1 

+ - m th (i) ujf h x'(i) I . (6) 



The coupling between nucleons and thermometer particles, H , is assumed 
to be weak, repulsive and of short range. It has to be as weak as possible 
in order not to influence the nuclear system too much. On the other hand 
it has to be strong enough to allow for reasonable equilibration times. Our 
choice is to put more emphasis on small correlation energies, smaller than the 
excitation energy, and to tolerate long equilibration times. 

The determination of the caloric curve is done in the following way. The nu- 
cleus is excited by displacing all wave packets from their ground-state posi- 
tions randomly. Both, centre of mass momentum and total angular momentum 
are kept fixed at zero. To allow a first equilibration between the wave-packets 
of the nucleus and those of the thermometer the system is evolved over a long 
time (10000 fm/c). After that a time-averaging of the energy of the nucle- 
onic system as well as of the thermometer is performed over a time interval 
of 10000 fm/c. During this time interval the mean of the nucleonic excitation 
energy 



-i N steps 

E* = ^— E ( Qniti) I H n I Qn(U) ) - E (N, Z) (7) 

^ steps i = i 

is evaluated. E (N, Z) denotes the FMD ground-state energy of the isotope 
under consideration. The time-averaged energy of the thermometer E t h, which 
is calculated during the same time interval, determines the temperature T 
through the relation for an ideal gas of distinguishable particles in a common 
harmonic oscillator potential (Boltzmann statistics) 
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T = four- 



th 
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'E th /N th + \huj th 
E th /N th - Ihuth, 



(8) 



3 The caloric curve 



The relation between the excitation energy and the temperature is determined 
for the four nuclei 16 0, 24 Mg, 27 Al and 40 Ca using the same container potential 
with ftiu — 1 MeV. In addition the dependence on oo is investigated for 16 
leading to an estimate of the critical temperature. 




E*/A=3.2MeV 
E*/A=7.9MeV 
EVA=10.8MeV 
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Fig. 1. L.h.s.: caloric curve of 24 Mg, 27 A1 and 40 Ca at ftiij = 1 MeV, r.h.s.: 
time-averaged radial density distri 
energies in the coexistence region. 



time-averaged radial density distribution of 24 Mg at various excitation 



The caloric curve shown in the graph on the left hand side of fig. 1 clearly 
exhibits three different parts. Beginning at small excitation energies the tem- 
perature rises steeply with increasing energy as expected for the shell model. 
The nucleons remain bound in the excited nucleus which behaves like a drop 
of liquid. At an excitation energy of 3 MeV per nucleon the curve flattens and 
stays almost constant up to about 11 MeV. This plateau at T ps 5 MeV has its 
origin in the coexistence of liquid and vapour phases, the latter consisting of 
evaporated nucleons which are in equilibrium with the residual liquid drop due 
to the containment. At the beginning of the plateau the system contains only 
few evaporated nucleons, towards the end more and also small intermittent 
condensed fragments which may amalgamate or dissolve into vapour again. 
Around E*/A m 11 MeV all nucleons are unbound and the system has reached 
the vapour phase. This is reflected by the rise of the caloric curve beyond this 
point. The "error bars" denote the r.m.s. fluctuations in temperature and ex- 
citation energy, respectively, which arise from the energy exchange between 
thermometer and nucleons. Only the total energy of the system (nucleons & 
thermometer) is conserved. The magnitude of the temperature fluctuations is 
larger than those of E* / A because the heat capacity of the thermometer is 
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smaller than that of the nucleus. 



One has to keep in mind that the plateau, which due to finite size effects 
is rounded, is not the result of a Maxwell construction as in nuclear mat- 
ter calculations. In the excitation energy range between 3 and 11 MeV per 
particle an increasing number of nucleons is found in the vapour phase out- 
side the liquid phase. This can be seen in the density plot on the right hand 
side of figure 1, where the radial dependence of the time-averaged density of 
a system of 24 nucleons is shown at three excitation energies in the coexis- 
tence region. For small excitations (E*/A = 3.2 MeV) the nucleus, which is 
bouncing around due to recoil, is surrounded with very low density vapour 
(solid line). The dashed-dotted line (E*/A = 7.9 MeV) and the dashed line 
(E*/A = 10.8 MeV) show that with increasing energy the vapour contribution 
is growing and the amount of liquid decreasing. However, in the high energy 
part of the plateau the averaged one-body density displayed here does not 
represent the physical situation adequately. The time-dependent many-body 
state shows the formation and disintegration of several small drops. Above 
E*/A 13 MeV only vapour is observed. 
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• 197 Au+ 197 Au, 600 AMeV 

□ 12 C, 18 + na, Ag, 197 Au, 30-84 AMeV 



▲ Ne+ Ta, 8 AMeV 



f (<E„>/<A > - 2 MeV) 



The caloric curve shown in fig. 1 has a 
striking similarity with the caloric curve 
determined by the ALADIN group [6] which 
is displayed in fig. 2. The position and 
the extension of the plateau agree with 
the FMD calculation using a containing 
oscillator potential of hoj — 1 MeV. Nev- 
ertheless, there are important differences. 
The measurement addresses an expand- 
ing non-equilibrium system, but the cal- 
culation deals with a contained equilib- 
rium system. In addition the used ther- 
mometers differ; the experiment employs 
an isotope thermometer based on chem- 
ical equilibrium and the calculation uses 
an ideal gas thermometer. One explana- 
tion why the thermodynamic description 
of the experimental situation works and 
compares nicely to the equilibrium result might be, that the excited specta- 
tor matter equilibrates faster into the coexistence region [16] than it expands 
and cools. The assumption of such a transient equilibrium situation [10,17,18] 
seems to work rather well at least in the plateau region. Further investigations 
on the FMD time evolution of excited nuclei without container will focus on 
such assumptions. 
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Fig. 2. Caloric Curve determined by 
the ALADIN group [6]. 



The thermodynamic properties of the nucleonic system are controlled by the 
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external parameter u similar to the volume in the case of an ideal gas. An 
increasing u leads to smaller volume and to higher pressure in the container 
which shifts the plateau of the caloric curve to higher temperatures and de- 
creases its extension, i.e. the latent heat. The critical temperature T c is reached 
for the uj c at which the plateau vanishes. Figure 3 displays this dependence 
for 16 O. The caloric curve is evaluated for three different oscillator frequencies 
of the container potential. A pronounced plateau is seen in the plot on the 
left hand side, where the oscillator does not influence the self bound nucleus 
very much. In the middle part the more narrow container potential is already 
squeezing the ground state, its energy goes up to E/A pa —5 MeV. The plateau 
is shifted to T pa 7 MeV and the latent heat is decreased. On the right hand 
side, for huj = 18 MeV, the coexistence region has almost vanished. In addition 
one observes very large fluctuations of T and E. 




E/A (MeV) E/A (MeV) E/A (MeV) 



Fig. 3. Caloric curve of O for the frequencies Hu = 1, 6, 18 MeV of the 
container potential. The solid lines show the low temperature behaviour 
of an ideal gas of 16 fermions in a common harmonic oscillator with level 
spacing huJ e Q , the dashed lines denote their high temperature behaviour 
in the confining oscillator (huj). 

The critical temperature T c , which can be estimated from the disappearance 
of the coexistence phase in figure 3, is about 10 MeV and has to be compared 
to the results of ref. [1,3,19] for finite nuclei including Coulomb and surface 
effects. All authors report a week dependence of the critical temperature on 
the mass number in the region from calcium to lead. Jaqaman's result with the 
Skyrme ZR3 interaction [1] can be extrapolated to 16 O to give T c pa 8 MeV, 
Bonche [3] arrives at the same number using the SKM interaction, but gets 
T c pa 11 MeV with the SIII interaction. Close to the last result is the value 
extrapolated from ref. [19] where T c pa 11.5 MeV for Gogny's Dl interaction. 
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